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Abstract. Neutral atoms may be trapped via the interaction of their magnetic dipole moment with mag- 
netic field gradients. One of the possible schemes is the cloverleaf trap. It is often desirable to have at 
hand a fast and precise technique for measuring the magnetic field distribution. We use for instantaneous 
imaging the equipotential lines of the magnetic field a diagnostic tool which is based on spatially resolved 
observation of the fluorescence emitted by a hot beam of sodium atoms crossing a thin slice of resonant 
laser light within the magnetic field region to be investigated. The inhomogeneous magnetic field spatially 
modulates the resonance condition between the Zeeman-shifted hyperfine sublevels and the laser light and 
therefore the amount of scattered photons. We apply this technique for mapping the field of our cloverleaf 
trap in three dimensions under various conditions. 

PACS. 32.60.+i Zeeman and Stark effects - 07.55.Ge Magnetometers for magnetic field measurements 



1 Introduction 



For two decades, the invention of cooling and trapping 
techniques for neutral atoms boosts the field of cold atomic 
physics and several years ago even led to the realization 
of the long expected goal of Bose-Einstein condensation 
(BEC) in dilute atomic gases, followed by a wealth of re- 
markable experiments on atom lasing and superfluidity of 
condensed gases. Most trapping configurations are based 
on light forces or magnetic field gradients or combina- 
tions of both. The most frequently used trap for neutral 
atom trapping is the Magneto- Optical Trap (MOT) 
MOTs have allowed the confinement of more than 10™ 
atoms |2j at temperatures well below 100 fjK and den- 
sities up to 10 12 cm" 3 . Radiation trapping by the opti- 
cally thick cloud, however, sets a limit to further com- 
pression due to photon rescattering ||. Even employing 
more sophisticated schemes j|, the highest phase space 
density achieved today remains at about one order of mag- 
nitude below the threshold to Bose-Einstein condensation. 
The phenomenon of radiation trapping is intrinsically con- 
nected to the fact that in magneto-optical traps the restor- 
ing force comes from the radiation pressure exerted by the 
laser beams. Alternative trapping schemes based on the 
dipole force of far-detuned laser beams avoid the radiation 
trapping problem |^,^|. A different approach is to exploit 
the weak interaction of the dipole moment of paramag- 
netic atoms with magnetic field gradients. Magnetic fields 



can be designed to generate local minima in free space 
capable of acting as trapping potentials. Several types of 
magnetic traps have been successfully used in BEC ex- 
periments like the Time- Orbiting Potential (TOP) trap 
M and various Ioffe-Pritchard (IP) type traps S, ||. Be- 
cause magnetic traps are generally weak, large magnetic 
fields and loading with precooled atoms is required. Gen- 
erally the atoms are precooled in a standard MOT and 
then transferred into the magnetic trap. In the absence of 
laser light the magnetic trap requires a new kind of cool- 
ing mechanism. The only mechanism that has successfully 
led to BEC up to now is Forced Radiofrequency Evapora- 
tion jlO],[ll]] . A necessary condition for evaporation to work 
is that the rate of elastic collisions between the trapped 
atoms be large enough to maintain the cloud in thermal 
equilibrium during the evaporation process which must 
take place within a period shorter than the trap lifetime 
[ p"2| . Large collision rates mean big atomic clouds and/or 
steep trapping potentials. 

This work is motivated by our need to characterize 
our cloverleaf trap which is, in principle, a Ioffe-Pritchard 
type trap |L3| . The necessity of optimizing the efficiency 
for loading, trapping and evaporative cooling sets tight 
conditions for the design of the magnetic coils and the 
current supplying circuitry. E.g., the field gradient must 
be rather strong and capable of being varied over a wide 
range. It is important to calibrate the secular frequencies, 
which essentially govern the shape of the trapped atomic 
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cloud, as a function of the control parameters. Therefore 
we need a precise control of the fields and a suitable way 
of measuring and monitoring them. Measuring magnetic 
fields is usually performed with Hall probes. However, very 
often the region of space to be analyzed is within a vac- 
uum recipient and thus not accessible to massive probes. 
We use the idea of fluorescence imaging magnetic fields 
jnj to develop a simple technique for the instantaneous 
mapping of magnetic field cquipotcntial lines in two di- 
mensions within a closed vacuum recipient. This technique 
that we call Magnetic Field Tomography (MFT) can even 
be extended to allow a three-dimensional reconstruction 
of the complete magnetic vector field. The method pro- 
vides a simple and reliable diagnostic tool for quantifying 
the field distribution of our magnetic trap for any set of 
values of the control parameters. 

This paper is organized as follows. In section 2, we 
present our magnetic trapping experiment and briefly ex- 
pose our approach to calculating the magnetic field distri- 
bution. In section 3, we introduce our Magnetic Field To- 
mography technique. We tested the technique and present 
experimental results in section 4, where we also show how 
to process the experimental data and compare the results 
to the calculated fields. We conclude our paper with a 
general discussion. 



2 The Magnetic Trap 

Magnetic traps are usually built with current carrying 
coils or permanent magnets [ fl5| . In the first case, high 
magnetic field gradients generally require high currents 
flowing through the coils. In special cases , hig h gradients 
may be achieved with lower currents [|lq,|l7|| . The main 
constraints for the design of the coils for the cloverleaf 
trap arise from the geometry of the vacuum chamber and 
by the power supply available: in our setup, the minimum 
distance between the coils and the center of the trap is 
15 mm. We apply up to 300 A and need to dissipate the 
10 kW resistive heat by cooling the coils with a high- 
pressure water flow passing through the hollow wires. Our 
magnetic cloverleaf trap is schematically shown in Fig- 
ure 1. Radial confinement of the atoms is assured by four 
pairs of coils in anti-Helmholtz configuration (off-center 
coils in Fig. 1) producing a quadrupolar waveguide field 
along the symmetry axis. Axial confinement is realized 
via a magnetic bottle (inner and outer axially centered 
coils in Fig. 1). The (inner) pinch coils create the neces- 
sary field curvature, and a pair of (outer) antibias coils 
in Helmholtz configuration compensates the large central 
offset produced by the pinch coils to a variable amount. 
As said earlier, to initiate run-away evaporative cooling, 
we must be able to realize large interatomic elastic colli- 
sion rates. This is achieved by compressing the magnetic 
trap: at the time of loading the trap from the MOT the 
mean secular frequency is typically lOtrap ~ 2-7T x 20 Hz 
(the transfer efficiency is best when the curvatures of the 
MOT and the magnetic trap are approximately matched), 
and the secular frequency of the fully compressed trap is 



Fig. 1. Scheme of the cloverleaf trap. The cloverleaf coils (blue) 
produce a radially linear magnetic field. The pinch coils (red) 
produce an axially confining magnetic bottle. The antibias coils 
(green) are Helmholtz coils and serve to compensate the mag- 
netic field offset. 



Utrap ~ 27r x 200 Hz. In the case of the cloverleaf trap, the 
compression is controlled via the antibias field. 

Precise control of the magnetic field bias is also impor- 
tant because it fixes an offset for the resonant radiofre- 
quency used for evaporative cooling: The evaporating sur- 
face |Q is set by the condition that the radiofrequency 
balances the shift of the atomic energy levels induced by 
the local magnetic field T%Lo r f = |/xB(r)|. In practice, the 
magnetic trap is compressed until the bias amounts to 
only a few Gauss. This means that fluctuations should 
be very small which is not trivial because the bias results 
from the subtraction of two large magnetic fields, the mag- 
netic bottle field and the antibias field. For our setup this 
calls for a current stability superior to 10 -4 . Furthermore, 
when loading the magnetic trap or releasing the atoms, we 
need to be able to quickly switch on and off the magnetic 
field, i.e. typically faster than within 100 /is |ll[. In sum- 
mary, the current control circuitry has to fulfill several 
demanding requirements, since we ask for high current, 
low noise, independent tuning of part of the current, and 
fast switching. The details of the current switch will be 
detailed elsewhere 

Once the geometry of the coils is known, calculating 
the magnetic fields B(r) is very easy. For simple geomet- 
rical shapes of the coils, one can use analytic formulae 
[ p0[ containing elliptical integrals. For complicated shapes 
of the coils or asymmetric arrangements, we numerically 
integrate the Biot-Savart formula 



B(r) 



fiol f ds x (r — s) 
|r - s| 3 ; 



(1) 



where s = s(<r) is the current path parametrized in small 
steps <f. From such calculations, we determine to first order 
the following data characterizing our cloverleaf trap close 
to its center. The magnetic field gradient per ampere pro- 
duced by the current I c i ov flowing through the cloverleaf 
coils is d r B c i ov /I c i ov = 0.438 G/cm/A. The curvature of 
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Fig. 3. Level scheme of the sodium D2-line between the hy- 
perfine levels F = 2 and F' — 3. The axial magnetic field B z 
is assumed linear in the coordinate z. 



Fig. 2. Magnetic trapping potentials in the isotropic and the 
compressed case. The experimental parameters only differ in 
the amount of antibias current: I an ti — 50 A in the isotropic 
case and Ianti — 85 A in the compressed case. The planes slic- 
ing the potentials represent the laser being tuned to a Zeeman 
shifted resonance. Our MFT technique images the contour lines 
are at these detunings (see text). 



the magnetic bottle field generated by the current I p i nc h 
in the pinch coils is in radial direction B P i nc h / Ipinch = 
—0.186 G/cm 2 /A and in axial direction d^B p i nc h/ Ipinch — 
0.373 G/crn/A. And the magnetic field amplitude pro- 
duced by the pinch and the antibias fields at the trap 
center are B pinch /I pm ch = 0.767 G/A and B anti /I anti = 
2.693 G/A. We usually work with the currents I c i ov = 
Ipinch — 285 A and I an ti variable. Inaccuracies in the geo- 
metric shape of the coils limit the level of precision of the 
calculations to 10% uncertainty. However, as we will show 
later, we can calibrate the calculations using the MFT 
method and achieve uncertainties lower than 2%. 

The absolute value of the magnetic fields |B(r)| is 
shown in Figure 2 for two different values of the antibias 
current. For low antibias current (left side), I an ti = 50 A, 
the trap is roughly isotropic at the center, the mean cur- 
vature of the magnetic field is small, and there is a large 
magnetic field offset, B(0) — 100 G. However, when we 
apply a high antibias current, I an ti — 85 A, the trap is 
axially elongated, radially compressed, and the magnetic 
field offset is small, B(0) = 1 G. 



3 Magnetic Field Tomography 

The basic idea of magnetic field tomography can be sub- 
sumed as "laser excitation spectroscopy on Zeeman-shifted 
electronic transitions with spatially resolved fluorescence 
detection" . We are detailing this idea in the following. The 
Hamiltonian of atoms with the electron total angular mo- 
mentum jij = gj^B^ and the nuclear spin fij = gi^i^l 
subject to a magnetic field B reads pif 



Hhfs + Hb 



6(IJ) 2 + 3IJ-2/(/+l)J(J+l) 
_ A,U + Bj 2/(21 -1)2J(2J-1) 

-/1,/B - /x/B 

The first term corresponds to the magnetic dipole in- 
teraction, and Aj is called the hyperfine constant. For 
sodium it takes the values Aj(3 2 Si/ 2 ) = 885.8 MHz 
and Aj(3 2 F 3/2 ) = 18.65 MHz [|l|. The second term 
corresponds to the electrostatic interaction of the nuclear 
quadrupole moment. The constant of the quadrupole in- 
teraction for sodium, Bj(3 2 P^/2) — 2.82 MHz, is small 
but not negligible. The next two terms describe the inter- 
action of the atom with the external magnetic field. The 
contribution containing the nuclear magneton /Xjv is neg- 
ligibly small. It is interesting to point out the behavior of 
the fully stretched spin states \F = J + 1 , nip = ±F). For 
those states the shift due to the interaction (Q) is linear 
in |B|. The shifts of all other states, since our magnetic 
trap operates in an intermediate regime of magnetic fields 
where the hyperfine coupling is perturbed by the mag- 
netic fields, are conveniently numerically calculated as the 
eigenvalues of the Hamiltonian. In order to keep the anal- 
ysis simple, we will first concentrate on the fully stretched 
states. However, we will need to stress the full expression 
to describe several features of our observations. 

When a laser is irradiated on an electronic transition 
with Zeeman degeneracy, the situation gets more compli- 
cated, because the Zeeman substates can be coupled by 
the laser, optical pumping takes place, and the energy 
levels are additionally shifted by the dynamical Stark- 
effect by an amount that depends on the Rabi frequency 
fl = dE/h, where d denotes the electric dipole moment of 
the transition, and E denotes the electrical field amplitude 
of the light. Wc will focus here on an approximation that 
holds in our experiment, i.e. the Rabi frequency is weak 
compared to the Larmor frequency top = /xB/7i. Then the 
energy levels are mainly shifted by the Zeeman-effect, and 
we may disregard the light shift. In this case, the magnetic 
field direction lends itself naturally as the quantization 
axis, because the Hamiltonian is diagonal with respect to 
this axis. 

For a laser irradiated on an electronic resonance to in- 
duce transitions, these transitions must be allowed by se- 
lection rules. Electric dipole transitions e.g. require F' — 
F = 0, ±1 and m F > — mp = 0, ±1, where the prime 
denotes the excited level. In an external magnetic field 
these transitions are not degenerate, so that magnetic field 



4 



Ph.W. Courteille et al.: Magnetic Field Tomography 



probing via observation of the fluorescence as a function 
of laser detuning is not unambiguous. If we choose our 
atomic transition so that the total angular momentum of 
the upper level is superior by 1 to that of the lower level, 
F' = F + 1, and if we apply a particular laser polar- 
ization, i.e. a + or a~ respectively, we drive a so-called 
cycling transition, i.e. the excited state can only decay 
into the very same ground state (see Fig. 3). The cycling 
transition involves the levels Si/ 2 ,F — 2,mp — ±2 and 
P 3 / 2 , F' — 3, mpi = ±3. Assume for example a laser beam 
irradiated along the z-axis, k = ke z , inside a magnetic 
field oriented in the same direction, B = B z e z . The choice 
of the z-axis defines the labelling of the light polarization, 
i.e. the light polarization called a~ is the one that excites 
transitions with mpi — mp = —1. Let us consider a o~ 
polarized red-detuned laser A < 0. For atoms subject to 
a positive magnetic field, B z > the laser then drives the 
cycling transition mp = F = 2 and mpi — F' — 3, while 
atoms that are subject to negative B z < are off reso- 
nance. By inverting either the polarization or the detun- 
ing, we probe the atoms that are subject to negative B z . 
This simple picture only holds if we can assume to have 
a strong magnetic offset field in z-direction, |B| « \B Z \, 
which is the case for the field to be investigated here. We 
will see that the condition is less restrictive than it seems 
in the first place. 

Because our cycling transition effectively realizes a two- 
level atom, the expression for the resonance fluorescence 
in terms of scattered photon numbers S per probing time 
tprobe takes a very simple form: 



S — — o~N a ttp ro i, e , 
nui 



(3) 



where I is the intensity of the irradiated light, N at the 
number of atoms in the interaction zone, and the optical 
cross section is 



a r 2 



4[/\-7r 1 (/x'-^) B ]' 



2Q 2 + r 2 



(4) 



where r = 2ir x 9.89 MHz is th e natural linewidth of 
the sodium D2 line, and Q = \J o^FI /huj is the Rabi 
frequency. For the cycling transition between the fully 
stretched states we find (^'—^)B = /is(g_F'm_F' —gpmp)s\Ti 
±/Lts<r|B| for (T~^~ polarized light, where = sign (B z ) and 
the Lande factors are gp = 1/2 and gpi — 2/3. 

The light is scattered at a nondegenerate transition 
between magnetic sublevels so that we have to weight the 
coupling strength {i.e. the Rabi frequency) with the rela- 
tive strength of the specific transition. The relative cou- 
pling strength can be expressed by {6j} symbols which 
reflect the fine and the hyperfine structure coupling and 
by (3j) symbols for the coupling of the atomic angular 
momenta to the magnetic field |2lJ . For the cycling tran- 
sition of interest, the relative coupling strength is 1. In 
those cases where B is not perfectly parallel to e z , the 
coupling strength depends on the relative orientation of 
the polarization of the irradiated light and the magnetic 
field direction eB. Those effects can be included in our de- 



Fig. 4. Scheme of the tomographic measurement experiment 
with atomic beam, laser light sheet and imaging system. The 
symmetry axis of the cloverleaf trap is along e z . The imaging 
system is inclined by 22.5° with respect to the atomic beam 



scription by using a modified Rabi frequency How- 
ever, we will see later that we do not need to specify the 
functional dependence for our purpose. 

The above considerations hold for single atoms inter- 
acting with a magnetic and an optical field. In the fol- 
lowing, we will address the question how these atoms can 
be utilized for magnetic field tomography. We consider 
spatially inhomogeneous fields B(r) and E(r) and specify 
the physical situation further by assuming that the spa- 
tial extend of the magnetic field is much larger than the 
interaction zone of the atoms with the light field. In fact, 
we will use a hot atomic beam perpendicularly intersected 
by a thin sheet of laser light. This allows us to disregard 
the impact of the magnetic field gradient on the atomic 
motion. We define our coordinate system by calling the 
propagation direction of the laser beam z-axis and the 
normal vector of the light sheet plane y-axis (see Fig. 4). 
The light sheet may be displaced along the y-axis by an 
amount yo and has an inhomogeneous Gaussian intensity 
distribution along the x-axis so that the intensity can be 
written 

I(x,y) = 2P/ttw 2 e~ 2x2/ < Xly0-d/2,y0+d/2](y), (5) 

where wo w 1 cm is the beam waist and P « 100 mW 
the laser power. The laser light has a well-defined polar- 
ization or 7T. The atomic beam is assumed to provide 
a homogeneous atomic density in the interaction region 
n{r) ss n ss 9 x 10 6 cm~ 3 . 

The fluorescence is collected along the y-axis, i.e. we 
integrate the fluorescence signal across the thin slice. A 
lens system images the scattered light so that we obtain a 
spatially modulated photon flux which can be expressed 
through a generalization of equation (|^) by 



Sa,vo,<;( x i z ) 

f d/2 I(x,y) 
at probe I ^ n(r)er(r, A)dy 



(6) 



-d/2 



Too 

at probe n drf2 S £.(x,yo,z) 2 



4[A-h-\[i B \B(x,y 0) z)\] +2f2 iti (x,yo,z) 2 + r 2 

The z dependence of the Rabi frequency comes from the 
dependence of the coupling strength on the relative orien- 
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tation of £ and B. The photon flux arriving at the image 
plane is recorded with a CCD camera. The proportional- 
ity constant a can be estimated knowing the camera pixel 
size dp X i = 9 /zm, its sensitivity (including quantum effi- 
ciency and electronic gain) tjccd = 1 count/100 photons, 
the light collection solid angle of the imaging lens system 



f2 sa = 0.06%, and its magnification 



0.17, 



d 2 

u pxl 



VCCD^sa- 



(7) 



We typically integrate the signal for t pro b e = 500 ms. 

Equation (^|) immediately reveals that, because we as- 
sumed strong magnetic fields, the fluorescence nearly van- 
ishes everywhere where the resonance condition A = ±^ _1 /is|B(r) 
is not satisfied, or in other words, the laser light slices the 
magnetic field at a given field strength |B(r)| = ^phA/fiB- 
The spatial distribution of the fluorescence reflects the 
modulus of the magnetic field. This situation is depicted 
in Figure 2 by a plane that interpenetrates the magnetic 
potential. For a given detuning A, light sheet position 
yo j and polarization A we expect the fluorescence image 
S A . vo-s( x ' z ) to form a narrow ridge along a closed path 
which basically traces a two-dimensional equipotential con- 
tour line of the magnetic field constraint to the plane of 
the laser light sheet (see Fig. 5). At the top of the fluores- 
cence ridge the photon count rate reads 



n a drf2 i& (x,y ,z) 2 
2f2 iti (x,y ,z) 2 + r 2 ' 



(8) 



where TiA — shbB(t). We see that the location of the 
ridge is not influenced by the dependence of the Rabi fre- 
quency on eB, but the height of the ridge is modulated. 
Furthermore, the contrast of this modulation vanishes as 
saturation is approached, J? — > r. Another reason for a 
modulated fluorescence is the imperfect homogeneity of 
the atomic density n(r). This contribution, however, is 
easily avoided by normalizing the fluorescence image with 
an image at zero magnetic field. 

As we stated earlier, the cloverleaf trap is, in axial 
direction, a magnetic bottle, i.e. the trapping potential 
is situated on top of a large magnetic field offset point- 
ing in z-direction. In this case, £? Z /|B| « +1, and we get 
the whole equipotential line using a single a + (or o~) 
laser polarization. However, this does not hold true for 
very compressed traps when the field offset is compen- 
sated by a homogeneous antibias field. Then in some re- 
gions of space, we may have £? 2 /|B| < 0. For the same 
reason, quadrupolar fields (or more generally, fields which 
are antisymmetric with respect to the z-axis) are only half 
imaged with a single circularly polarized laser. However, 
if we irradiate linearly polarized light, which is a linear 
superposition of left and right circularly polarized light, 
we get the complete image ±1 Sa,vos{ x > z )- 



4 The Experiment 

In our experiment, we use a fast atomic beam to probe 
the magnetic field. Hot sodium atoms are ejected from a 



Fig. 5. (a) Calculated and measured fluorescence image of a 
light scattering atomic beam inside a cloverleaf-shaped mag- 
netic field with Id ov — 285 A and Ianti = 80 A. The laser 
is linearly polarized and detuned +50 MHz from resonance, 
(b-d) Equipotential lines derived from fluorescence images as 
shown in Fig. (a). The total current for the magnetic trap is 
Iciov = 285 A. In Fig. (b) the antibias current is I an u = 85 A, 
the detuning is +50 MHz, and the light is a~ polarized. The 
different lines correspond to different locations of the laser slice 
displaced from one another by 1 mm. In Fig. (c) the slice is 
kept at the center, the antibias current at Ianti ~ 80 A, and the 
detuning is ramped from +10 MHz to +50 MHz. In Fig. (d) the 
slice is in the center, the detuning is fixed at +25 MHz, and the 
antibias current is ramped from Ianti = 80 A to Ianti = 100 A. 



heated oven through a fine nozzle. After a 2 m long dis- 
tance, they arrive at the interaction region as a 2 cm wide 
highly collimated and homogeneous beam: We verified 
that the radial velocity distribution v±_ is narrow enough 
to allow us to neglect the transversal Doppler effect at the 
D2 atomic resonance, kv± < r. The atoms enter the mag- 
netic field with a velocity satisfying dtB/B <C Ql which 
allows the assumption that the atomic spin adiabatically 
follows the fields. 

In the interaction region the atoms cross a beam of 
near-resonant laser light. When magnetic fields are present, 
the fluorescence scattering is spatially modulated by the 
inhomogeneous Zeeman-shift of the resonant levels. The 
fluorescence can conveniently be imaged to a CCD cam- 
era. Blurring of the images by column integration along 
the imaging direction can be avoided by only irradiat- 
ing a thin sheet of light whose plane is orthogonal to the 
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Fig. 6. (a) Multiple magnetic field tomographic lines. When 
the dynamic range of the magnetic field strength is large 
enough, we observe the Zeeman splitting of the a + transitions 
as multiple tomographic lines. The figures correspond to the 
cloverleaf trap parameters I c i ov = 285 A and Ianti = 60 A, 
and the detuning is +20 MHz. The four ridges correspond from 
the inside to the outside to the transitions (mF,m' F ) — (2,3), 
(1,2), (0, 1), and ( — 1,0). The noise in the measured date has 
been reduced to increase the visibility of the features, (b) and 
(c) The fluorescence light is polarized, and the polarization de- 
pends on the local magnetic field. In the figures (b) it has been 
filtered with a polarizer oriented to transmit light being lin- 
early polarized along e = e z . In the figures (c) i = e x light is 
left through. 

imaging axis. In this way, two spatial dimensions of the 
magnetic field are directly imaged and the third is probed 
by moving the light sheet back and forth. The strength of 
the magnetic field is measured by tuning the laser light to 
the Zeeman-shifted resonances, thus slicing the magnetic 
potential at a depth that can precisely be set by the laser 
frequency (see Fig. 2). 

Our experiment is sketched in Figure 4, and an ex- 
ample of a fluorescence image in shown in Figure 5(a). 
For this image we operated our cloverleaf trap at full cur- 
rent I c i ov — 285 A, applied I a — 80 A antibias current, 
and located the laser light sheet at the center of the trap. 
The ring-shaped border corresponds to atoms being res- 
onant with laser light blue-shifted by A = 2ir x 50 MHz. 
Therefore, the magnetic field along this border line is 
|B(r)| = HA/hb = 36 G. 

By recording images of the equipotential lines for var- 
ious detunings and various positions of the laser light 
slices along the imaging axis, we are able to reconstruct 
the magnetic field in three dimensions. To demonstrate 
this, we record several images like the one shown in Fig- 
ure 5(a) with varied parameters. From the fluorescence 



maximum of every recorded image, we extract the closed 
ridge. For this purpose, we have written a data analyzing 
program which, from an arbitrary starting point, is capa- 
ble of climbing the closest maximum and discerning the 
path to follow the ridge. We plot all ridges into a common 
picture. For example, the right hand side of Figure 5(b) 
shows the equipotential lines for A — 2ir x 50 MHz and 
Ianti = 40 A at various positions of the slice, thus form- 
ing an equipotential surface, |B(r)| = const. If we know 
the surfaces for all values of const, we can reconstruct the 
three-dimensional magnetic vector field B(r). To record 
equipotential surfaces at other magnetic field strengths, 
we repeat the procedure for different laser detunings. This 
is demonstrated in Figure 5(c) where we kept the light 
sheet at the trap center and varied the detunings between 
10 MHz and 50 MHz. Figure 5(d) shows the equipotential 
lines for fixed detuning of A = 2ir x 50 MHz and centered 
light sheet position but varied antibias current. 

The left hand sides of the Figures 5(a-d) show the re- 
sult of calculations using formula (1). In order to improve 
the calculations, we calibrated them by measurements us- 
ing a procedure to be detailed in the following paragraph. 
The remaining small discrepancies between theory and ex- 
periment are probably due to the inhomogeneous intensity 
distribution within the laser beam. 

When we drive the D2 line with <r + light, five tran- 
sitions between the hyperfine levels 2 Si/2, F — 2 and 

2 -f3/2i F — 3 satisfy the selection rule mpi — uif = 1. 
Under certain conditions, we experimentally observe all 
of them simultaneously, since the laser can be resonant 
with different transitions simultaneously at different re- 
gions of space, where the magnetic field has different val- 
ues |B(r)| = TiA/ [Ats(<7F'TOF' — gFmp)]. The right hand 
side of Figure 6(a) shows such a fluorescence image. When 
performing the calculations (left hand side of the figure), 
it is essential to take into account the nonlinear depen- 
dence of the Zeeman-shift on the magnetic field for transi- 
tions between non-stretched states (Paschen-Back regime 
of the hyperfine structure) . For calculating the curves in 
Fig. 6(a), we numerically diagonalize the matrix (Q) for 
the ground and excited hyperfine level. 

Furthermore, we observe that the fluorescence light is 
polarized and that the polarization of is spatially mod- 
ulated by the magnetic field. For the images shown in 
Figure 6(b) and (c), we filtered the fluorescence light in 
front of the camera with a linear polarizer. The fluores- 
cence ridge does not move, however, its contrast is strongly 
modulated and changes when we set the transmitting axis 
of the polarizer along n = e z (Figure 6(b)) or n = e x (Fig- 
ure 6(c)). The behavior is easily understood in terms of 
the radiation pattern of an atomic dipole inside a magnetic 
field. In all directions orthogonal to B the atom emits lin- 
early polarized a light. In the plane which is sliced by our 
laser light sheet, our cloverleaf trap yields B(a;, 0, z) _L e^, 
so that we expect linearly polarized light, indeed. And 
since B(x, 0, z) rotates along the fluorescence ridge, the 
polarization of the scattered light does so, as well. The 
left two pictures show the result of calculations according 
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Fig. 7. Cuts through fluorescence images along a line defined 
by x — y = for various laser intensities I = 5.6, 12, 118, and 
186 mW. The effect of power-broadening is clearly visible. The 
magnetic field parameters were I c iov = 285 A, Ianti = 85 A 
and +40 MHz. 

to 

The resolution of the MFT method can be defined as 
the halfwidth of the fluorescence ridge observed at the 
CCD camera. With this definition, the width 2 Sz of the 
emerging contour line in e z direction immediately fol- 
lows from equation (|), 26 z = hJr 2 + 2Q^/((j, B d z \B\). 

Thus it depends on the magnetic field gradient and on 
the (power-broadened) linewidth of the D2 resonance. In 
terms of magnetic fields the resolution is given by SB = 
<9j|B| 5z, and we see that, even if we work below satura- 
tion, Q <C r, we are not sensitive to magnetic field vari- 
ations that are weaker than SB w 3.6 G. Figure 7 shows 
cuts through fluorescence images along a line defined by 
x = y = for various laser intensities. The fluorescence 
rate increases proportionally with low laser intensities, but 
saturates at high intensities. The highest spatial resolu- 
tion found was about 1 mm 3 at a location where the axial 
field gradient was c^B sa 50 G/cm which corresponds to 
SB w 5 G. We could, in principle, use narrower (Raman) 
transitions. However, the fluorescence rate will be lowered 
as well, thus requiring a more sensitive detection system. 

The MFT method can also be used to precisely deter- 
mine the amount of antibias current necessary for exact 
compensation of the offset field created by the pinch coils. 
For this purpose, we set the antibias current to some value 
Ianti, observe the fluorescence distribution (Fig. 5(a)) on- 
line, and find the laser detuning A at which the fluores- 
cence disappears. This gives us the minimum of B z . We 
repeat this procedure for various antibias currents and ob- 
tain two curves, one for o~ and one for a + polarized light. 



Fig. 8. Compensation of the magnetic field offset via I an ti 
using I c i ov = 285 A. Shown is the detuning at which atoms at 
the bottom of the potential scatter light. These data can be 
used to calibrate the calculations of the magnetic fields. 

The curves intersects near A = 0. At this detuning the an- 
tibias is perfectly compensated, i.e. the magnetic trap has 
maximum compression. The curves are shown in Figure 8. 
The compensation current is I an ti = (86.1 ± 1.1) A. As we 
see, the uncertainty of this method is rather small. Fur- 
thermore, the slopes of the curves in Figure 8 yield the 
calibration of the pinch coils, B pnc h/I c iov — 0.830 G/A, 
and of the antibias coils, B anti /I ant i = 2.743 G/A, at the 
center of the trap. We can now utilize these data to cor- 
rect our magnetic field calculations based on Formula ([!]). 
We also calibrated the current I c i ov using the same method 
and obtained values for the cloverleaf trap parameters that 
agree well with the above ones. 



5 Conclusion 

In conclusion, magnetic fields have been probed using a 
technique well-suited for fast three-dimensional tomographic 
reconstruction of a magnetic field within a volume which 
is not accessible to material detectors like Hall probes, 
e.g. within an ultra-high vacuum. The drawback of this 
method is that the resolution is obviously limited by the 
natural linewidth of the transition, so that it is only appli- 
cable for large magnetic fields and large field gradients. A 
solution to this problem could be the use of radiofrequency- 
optical double resonance schemes [p2|,p3|. 

The MFT method is well adapted to investigating mag- 
netic trapping fields for neutral atoms. In particular, it 
permits precise mapping of the trap geometry and fast lo- 
calisation of its minimum. The response time of the MFT 
method is only limited by the time required for optical 
pumping. Therefore the method could be extended for 
monitoring rapid changes in the magnetic field strength 
being as fast as a few 10 ns using stroboscopic imaging or 
even for real-time control of the magnetic trap configura- 
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tion. We conclude that, since this technique only requires 
a hot atomic beam and a weak tunable laser beam, we be- 
lieve it to be a simple and versatile diagnostic instrument 
for probing magnetic fields. 
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